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Abstract. A novel definition of the Laplacian of 1/r is presented, suitable for 
advanced undergraduates. 



1. Introduction 

Discussions of the Laplacian of 1/r generally start abruptly, in medias res, by stating 
the relation 

V 2 ^ = -An5 3 (r), (1) 

where r is the magnitude of radius vector r and S 3 (r) is the three-dimensional delta 
function, which is then proved in various ways, clarifying thus its meaning. A glance 
at equation (1) reveals, however, that the symbol V 2 appearing in it can not have its 
ordinary, classical meaning of V • V, where, in Cartesian coordinates, 
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is Hamilton's operator nabla (cf, e.g., Redzic [1]), since the classical expression V 2 (l/r) 
vanishes for r ^ and is not defined at r = 0. Therefore, instead of the familiar form 
(1), in this note we will use a less confusing notation [2J [3] 

V 2 - = -4vr5 3 (r); (2) 
r 

the expression V 2 (l/r) we will call the generalized (distributional) Laplacian of 1/r and 
try to fathom its meaning. Let us briefly review some typical proofs of (2). 

The well-known way to demonstrate (2) is to regularize 1/r in terms of a parameter 
a so that the regularized function is well-behaved everywhere for a ^ 0. Then verification 
of (2) consists in showing that in the limit a — > 0, — 1/Att times the Laplacian of the 
regularized function is a representation of the three-dimensional delta function 5 3 (r). 
For example, regularizing 1/r as l/\/r 2 + a 2 , Jackson [1] shows that 

V 2 - = lim V 2 ; 1 = -An5 3 (r); (3) 
r a^o Vr 2 + a 2 
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the limit here is the weak limit (cf, e.g., [5]). A more sophisticated method of proving 
(2) would be to derive first the generalized second-order partial derivatives of 1/r with 
respect to Cartesian coordinates [3j E 
A generalization of (2) 
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obtained simply by a different choice of the origin, can be demonstrated either via the 
regularization procedural or by employing a well-known electrostatic argument. Namely, 
it can be shown that the 'potential' 
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where the 'density' g{r') plays the role of a well-behaved 'test' function, satisfies 
Poisson's equation 
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= -4vrf?(r), (6) 

by making use of the divergence theorem and Gauss's theorem from electrostatics [8j [9] . 
Since the right-hand side of equation (6) can be written as 

-An J g(r')5 s (r-r')d 3 r', (7) 

it follows that the operator V 2 can enter an integral of the form (5) under proviso 
that 'during entrance' it converts into the generalized operator V 2 whose action on the 
function l/\r — r'\ is defined by equation (4). 

On the other hand, the representation formula for a well-behaved scalar function 
of position $ 
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f A typical informal derivation of (2) (cf, e.g., [7]) starts from relation 



/V*id»r = /V.(-^ 



d d r, 



which, using the divergence theorem, equals —An — —An J S 3 (r)d 3 r. However, the use of the divergence 
theorem is not legitimate here, since the function — r/r 3 is singular at r = 0, and the result 



V -d d r 
r 
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is incorrect since the left-hand side of the last equation equals zero. 
| Jackson [3], in fact, shows the validity of a generalization of (3) 
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where the point r is within the volume V and 1Z = \r — r'\, is obtained from Green's 
second identity 

where the volume V — V E is obtained by excluding a small ball of radius e and centre 
at r' — r from the volume V, and is the surface of the ball, taking the limit e — > 
(cf, e.g. [ID]). In this way, the singularity of the function 1/TZ at r' — r is managed 
with, making possible the use of Green's second identity. As is well known, this classical 
procedure provides grounds for replacing the term V' 2 {1/71) in (9) by its generalized 
counterpart V' 2 {\/1Z) defined by (4), extending of course at the same time the restricted 
(singularity-free) volume of integration V — V e to the whole region V and removing the 
integral over the surface S £ \^ Thus, employing of the generalized Laplacian of the 
function l/\r — r'\ instead of the classical one in regions that contain the singular point 
r' — r, replaces lengthy procedures based on classical analysis, providing a shortcut to 
the correct final results. 

In this note, an alternative definition of the generalized Laplacian of 1/r will be 
presented. Taking into account the ubiquity and importance of this somewhat tricky 
concept, the alternative derivation of relation (2) could perhaps be of some pedagogical 
interest. 
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2. Integral versus differential definitions of classical operators 

A perusal of the literature reveals that a common feature of discussions of the Laplacian 
of 1/r is that, right from the outset, the Laplacian is identified with V 2 = V • V. 
That is, in the primary definition of the Laplacian of a scalar field, which is the 
divergence of the gradient of the field, both the divergence and gradient operators 
are understood according to their differential definitions, V- and V, respectively, 
in accord with common practice of defining familiar operators in vector analysis by 
differential operations (cf, e.g., [Hill]). However, as is well known, there exists also an 
alternative way of defining the classical operators by means of integral operations (cf, 
e.g., [HI EES H3] ) . While the integral and differential definitions are equivalent in the case 
of differentiable fields, the former definition provides physical insight and computational 
convenience, and thus should be preferred, as Sommerfeld suggested [12]. Moreover, it 
appears that the integral definitions of classical operators are more to the point than 

§ To put it rigorously, in this case relation 

lim(/ $( r ')V' 2 ^dV+ / [A|i-$A(l^l d( s'\ = _ 47r $( r ) 
*->°Uv-v. ft JsAKdn' dn'\n)\ J K) 

serves as a definition of V 2 (l/7?.) (cf, e.g., 5J, since 

-47r$(r) = -4tt J $(r')£ 3 (r' - r )d 3 
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the differential ones when applied to a singular point of the field. Let us examine the 
last point in some detail. 

As a simple illustration, consider the Laplacian of 1/r at the singular point r = 0. 
Using the differential definitions, V 2 (l/r) is not defined at r = 0. On the other hand, 
using the integral definition of the divergence of a vector field A at a point P, 

divA = lim- I A-dS, (10) 



where the flux of A is through a closed surface surrounding P and r is the enclosed 
volume, and setting A = grad(l/r) = — r/r 3 , we obtain that 

divgrad(l/r) = — 47roo = — oo, (11) 

at r = 0. Since infinite value of a function at a point is not permitted in classical 
analysis, the two definitions seem to be in a dead heat. However, this is not so, as the 
following argument will show. 

Calculate the average value of the Laplacian of 1/r over the volume of a ball of 
radius e and centre at r = 0. Using the differential definition of the Laplacian we obtain 

<V 2 (l/r)> = j^j^ j y V 2 (l/r)d 3 r = 0, (12) 

where (...) stands for the average value, and V £ is the volume of the ball. Once again, 
the same result is obtained using the integral definition of the Laplacian, 

(divgrad(l/r)) = 1 / div(-r/r 3 )d 3 r = 0, (13) 

since according to both definitions the Laplacian of 1/r vanishes for r ^ 0, and the value 
of integral is not affected by a countable number of singular points of the integrand. 

Inspection of equation (10) reveals, however, that a reasonable physical definition 
of the volume average of the divergence of a vector field A over a volume r would be 

(divA) = - <t A ■ dS, (14) 
T Js T 

where S T is the surface of the volume r, rather than the standard definition 

(divA) = -^divAd 3 r. (15) 

Namely, while both definitions yield the same result when the field A is differentiable 
over the volume r, they give different results when A is singular at a point inside r. For 
example, if E is the electrostatic field of a point charge q located at the origin, and r is 
a volume containing the origin, (div.E) would be equal to q/reo according to definition 
(14). On the other hand, according to definition (15), (divE) would be equal to zero, 
since divi? is zero for r ^ and, classically, is not defined at r = 0. Keeping in mind 
that the divergence of a vector field is a measure of the strength of a source or sink 
of field lines, it is clear that the standard definition (15) would yield the absurd result 
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that the average strength of a point source over a volume containing the source is zero 
Thus, the integral definition (14) of the average divergence not only implies the integral 
definition (10) of the divergence at a point, but also is adequate for regions that contain 
singularities (point sources or sinks) of field linesQ 

3. A novel definition of the Laplacian of 1/r 

The above analysis provides the opportunity of introducing an alternative definition of 
the generalized Laplacian of 1/r. 

As a first step, define the e-Laplacian of 1/r, L e (l/r), in terms of a parameter e as 

T 1 J ([f (-r/r 3 ) ■ dS]/(4/3)7re 3 , if r<e 

L e - = div £ grad- = { [ ^\. 1 ' J/ Kl > (16) 
r r divgrad(l/rj = 0, if r > e. 

Thus, for r < e, L E (l/r) is the average divergence (as defined by the integral definition 
(14)) of the gradient of 1/r over the volume of a ball of radius e and centre at r = 0; 
for r > e, L e (l/r) is the classical Laplacian of 1/r. Definition (16) implies that 

L £ - = --Me-r) (17) 



r e 3 



where Q(x) is the Heaviside step function, 
6(x) = 



0, if x < 

1, if x > 0. 
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When £ — y 0, from equation (17) we obtain 

hW = { - 00 ' * r = ° (19) 
e^o £ r 1 0, if r^0. V ; 

On the other hand, integrating the product of L e (l/r) and a well-behaved 'test' function 
f(r), 

f(r)d 3 r = - [ -J(r)d 3 r, (20) 

Jr<e £ 

|| Note that for a given field, a source of field lines and a source of the field need not necessarily coincide. 
For example, the electric field lines of a point charge which is forever in uniform motion emanate from 
the present position of the charge, whereas the source of the field is the charge at the retarded positions 

(cf, e.g. rum]). 

Recall that a valid definition, inter alia, should be adequate (' definitio sit adequatd\ Recall also that 
there are other cases when the average value of a physical quantity can not be defined via an integral 
of the classical local values of the quantity. For example, the average velocity of a particle, is primarily 
given by (v) = At - / At, where Ar is a displacement of the particle during a time interval At, and not 
by (v) = (1/ At) J At v(t)dt. (The two definitions are equivalent only if r is a differentiable function of 
t.) Similarly, the average charge density over a volume AV^ is primarily (g) = AQ/ AV, where AQ is a 
charge inside AV, and not (g) = (f/AV) J AV g(r)d 3 r, if g(r) is the classical local charge density. Of 
course, if v(t) and g(r) are described by generalized functions (distributions), such as the Dirac delta 
function (which is the case when v(t) changes abruptly during At, or there are point charges inside the 
volume AV), no discrepancy arises between the two kinds of definitions. 
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using expansion of f(r) in a Taylor series around r = 0, and taking the limit e — > 
yields 



lim 

£->0 
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= -4tt/(0). (21) 
As is well known, result (21) can be expressed as 

L- = wlimL e - = -47r5 3 (r); (22) 

r e^o r 

where a more suitable notation L for the generalized Laplacian is now used instead of 
V 2 and wlim stands for the weak limit (cf, e.g., [21[S])Lj Equation (22) is tantamount 
to equation (3), providing another definition of the generalized Laplacian of l/rjj 

Discuss now briefly a closely related problem of defining analogously the generalized 
charge density for a point charge q located at the origin. Obviously, the corresponding 
e-charge density should be defined as 

g/(4/3)vre 3 , if r < e 
g(r) = 0, if r > e, 



Qe(r) 



(23) 



which can be recast into 



3g 
Aire 3 



0(e -r). 



(24) 
(25) 



Passing details, we give the final result 

g(r) = wlim^ e (r) = q5 3 (r). 

£->0 

where g(r) is the generalized charge density. Note that lim e ^o / Q £ (r)d 3 r = q, whereas 
/[lim e _j.o Qe{. r )]d 3 r = 0. Thus, the volume charge density of a point charge q located at 
r = is naturally described by g E (r) and its weak limit g(r), in perfect analogy with 
the volume flux density of the flux of grad(l/r), which is naturally described by L e (l/r) 
and its weak limit L(l/r). 

+ As was pointed out in Section 2, the standard practice of using automatically V 2 for the classical 
Laplacian may be misleading. 

* Recall that the 2D analogue of equation (4) reads 

V 2 ln|s-s'| = 2ir6 2 (s-s'), 

where s and s' are 2D radius vectors and S 2 (s — s') is the 2D delta function [4j|8l[T5], which setting 
s' = yields 

V 2 lns = 2nS 2 (s), 

which is the 2D analogue of equation (3). To prove the last relation, instead of regularizing Ins in 
terms of a parameter a as In \/s 2 + a 2 , the alternative 2D definition of the generalized Laplacian of In s, 

l\ns = 27n5 2 (s), 

can be introduced, following a 2D procedure analogous to the 3D procedure described above. 



On the Laplacian ofl/r 



7 



In the same way we can define the generalized divergence of the electric field E(r, t) 
of a point charge q that is moving on a specific trajectory r g (t). Starting from the 
corresponding e- divergence 

. L l E-dS]/U/3)7T6 3 , if £<e 
divi? = 0, if £ > £, 



we 



where £ = |r — r 9 (t)|, using Gauss's law in its integral form, §^ =£ E ■ dS = q/e 
obtain that E(r, t) must satisfy equation 

d^E(r,t) = qSir - r ^ t)) =IM (27) 

where 

div-Efr, t) = wlim div E £(r, t) (28) 

is the generalized divergence of E(r,t). This is consistent with a general sine qua non 
that when the sources of a field are idealized as point, line or surface distributions of 
charge and/or current, described by generalized functions, generalized operators must 
be employed in the field or potential equations instead of the classical onesr 



4. Conclusions 



We presented a novel definition of the generalized Laplacian of 1/r, avoiding 
regularization of 1/r [2111], 'electrostatic' procedure [S][H], or Green's second identity [S]. 
The definition is constructed employing the integral definition of the divergence instead 
of the differential one; thus, the usual notation V 2 is replaced by a less suggestive one 
L. It is shown that the Laplacian of 1/r can be naturally construed as the volume flux 
density of the flux of grad(l/r), in the same way as the volume charge density of a point 
charge located at the origin, introducing a reasonable generalized density. We believe 
that our analysis provides a simple and insightful alternative to the earlier discussions 
of the concept, clarifying thus its meaning. 
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